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The trinification model is an interesting extension of the Standard Model (SM) based on the 
gauge group SU(3)c x SU(3)l x SU(3)r. We study its low-energy phenomenology by constructing 
a low-energy effective field theory, thereby reducing the number of particles and free parameters that 
need to be studied. The resulting model predicts that several new scalar particles have masses in 
the 0(100 GeV) range. We study a few of the interesting phenomenological scenarios, such as the 
presence of a light fermiophobic scalar in addition to a SM-like Higgs, or a degenerate (twin) Higgs 
state at 126 GeV. We point out regions of the parameter space that lead to measurable deviations 
from SM predictions of the Higgs couplings. Hence the trinification model awaits crucial tests at 
the Large Hadron Collider in the coming years. 

PACS numbers: 12.10.Dm 


I. INTRODUCTION 

The discovery of the Higgs boson [1, 2] marks the estab¬ 
lishment of the Standard Model (SM) of particle physics 
as the model that correctly describes physics at experi¬ 
mentally available energies to date. All SM particles have 
been discovered, and the experimental data gathered at 
particle colliders match the predictions of the SM to good 
precision [3]. Yet, the SM is regarded to be an incomplete 
theory of nature: it lacks a dark matter candidate and is 
incompatible with the observation of non-zero neutrino 
masses. Also, the fermion masses and mixings are free 
parameters that display hierarchical patterns, and parity 
violation is introduced by hand. Therefore our quest to¬ 
wards a better theory of nature requires us to extend the 
SM. 

Grand Unified Theories (GUTs) [4-8] are interesting 
extensions of the SM in which the SM gauge group 
SU (3 )c x SU (2)l x U(1)y is embedded in a larger simple 
gauge group. The exceptional group Eq is an attractive 
example of a GUT group [6-8]. It is anomaly-free and 
left-right-symmetric (LR-symmetric), and as such it pro¬ 
vides an explanation for parity violation in the SM by 
spontaneous symmetry breaking. It appears in the com- 
pactification of string theories, which leads to either four- 
dimensional Eq gauge symmetry or one of E§ s maximal 
subgroups [9, 10]. One of these maximal subgroups is the 
‘trinification group’ G333 = SU{3)c x SU(3)l x SU(3)r. 
Models based on the trinification group have been stud¬ 
ied in several contexts [11-18]. 

In this work, ‘trinification model’ will refer to the setup 
described in refs. [7, 19-23]. The setup described there is 
interesting for several reasons: fermion masses and mix¬ 
ings of the SM can be reproduced using only a few pa¬ 
rameters, with a satisfactory fit for the solar neutrino 
mass difference and the neutrino mixing pattern. Also, 
a Standard-Model-like Higgs with a mass close to 126 


* hetzel@thphys.uni-heidelberg.de 
t b.stech@thphys.uni-Heidelberg.de 


GeV appears in a large region of parameter space of the 
model. Furthermore, it gives predictions for the matrix 
element of neutrinoless double-beta decay and the neu¬ 
trino masses, which allow the model to be tested with 
low-energy experiments. It also allows for various inter¬ 
esting phenomenological scenarios, such as the presence 
of a light fermiophobic Higgs in addition to the Standard- 
Model-likc Higgs, or even a degenerate Higgs state at 126 
GeV. 

In order to compare the trinification model with ex¬ 
periment, a study of the low-energy phenomenology is 
necessary. Due to the large number of scalars, a study 
of the full scalar mass matrix is challenging. However, 
several of the scalar fields will obtain very large masses 
when the trinification symmetry is broken, and thus can 
be integrated out from the theory. The result is an ef¬ 
fective field theory with the LR-symmetric gauge group 
SU(3)c x SU(2)l x SU(2)r x U(\)b-l, and fewer scalar- 
fields than in the trinification model. This model has 
the same low-energy properties as the trinification model, 
but is easier to study. We will refer to this model as the 
low-energy trinification (LET) model. 

Neither the LET model nor the trinification model re¬ 
solves the hierarchy problem. This problem is hidden in 
the vacuum expectation values (vevs) used in the model, 
which are presently not understood. In our treatment, 
all dimensionful parameters and masses are fully deter¬ 
mined by these vevs multiplied by dimensionless coupling 
constants. Since these vevs are momentum- and scale- 
independent (except for wave-function renormalization), 
their use as fixed parameters is justified. 

Left-right symmetric models based on the gauge group 
SU(3) C x SU(2) l x SU(2) R x U(1) B -l [24-26] have 
been studied extensively in the literature. Moreover, 
these models have many features in common with the 
two-Higgs-doublet model (2HDM) [27]. However, the 
LET model has properties that distinguish it from more 
general LR-symmetric models and the 2HDM, due to 
the trinification origin at high energy scales. A LR- 
symmetric model in the context of the trinification model 
has not been studied before to the best of our knowledge. 
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Therefore the LET model merits a study. 

In this work, we explore the low-energy phenomenol¬ 
ogy of the LET model. We briefly describe the trinifi- 
cation model in section II, and subsequently derive the 
properties of the LET model in section III. In order to 
aid our study of the LET model, we first introduce the 
Single-Bidoublet (SB) model, a simplified form of the 
LET model that has fewer scalar fields. The possible 
phenomenological scenarios for both the SB model and 
the LET model are worked out in section IV. We show 
how these models may be distinguished experimentally 
from the SM by studying the modifications of the SM 
Higgs couplings in section V. It turns out that the LET 
model allows for interesting phenomenological scenarios, 
such as a very light fermiophobic scalar with a mass in 
the GeV range or a degenerate scalar state at 126 GeV. 
We discuss these scenarios in more detail in sections VI 
and VII. Our conclusions are presented in section VIII. 


II. THE TRINIFICATION MODEL 


The Higgs sector of the trinification model contains 
two complex scalar fields H t , H 2 in the (1,3,3) repre¬ 
sentation of G 333 = SU(3)c x SU(3 )l x SU(3)r. We 
can use the SU(3)l x SU(3)r gauge symmetry to bring 
the vev of H 1 into diagonal form: 


(Hi) 


1 (vi 0 0 \ 

- 7 = 0 h 0 , (H 2 ) 

V 2 V 0 0 Ml/ 


1 (V 2 0 0 

— ° b 2 63 

\ 0 M M 2 


(!) 

Here we employed a matrix notation in which SU(3)l in¬ 
dices run vertically and SU(3)r indices run horizontally. 
All vev parameters are taken to be real in order to avoid 
tree-level CP- violation. The second field H 2 is necessary 
to break the left-right symmetry of G 333 : it cannot be 
made diagonal once (Hi) is taken to be diagonal. The 
off-diagonal parameters M, b 3 are taken to be unequal, 
and thus break the left-right symmetry. We assume the 
presence of large hierarchies among the vev parameters. 
The parameters M\,M 2 ~ 10 13 GeV are of the order 
of the scale where the Standard-Model gauge couplings 
g 1 and g 2 unify. The off-diagonal vev M is taken to be 
an intermediate scale of order 10 10 GeV, but could also 
be as low as a few TeV. The gauge couplings gh,R of 
SU(3)l,r are equal above this scale, whereas below M 
the left-right symmetry is broken. The other vev param¬ 
eters contribute to the W-boson mass and are therefore 
much smaller than Mi, M 2 , M. As such, the former are 
constrained by the relation v\ + v 2 + b\ + b 2 + 6 2 = v 2 = 
(246 GeV) 2 . The used scalar potential is renormalizable 
and all its parameters are taken to be real in order to 
avoid tree-level GP-violation. 


The fermions are grouped into the fundamental rep¬ 
resentation 27 of Eq. They are two-component left- 
handed Weyl spinors with respect to the Lorentz group. 
The fermion field decomposes into a lepton field L , a 


left-handed quark field Qr, and a right-handed quark 
field Qr, which are assigned to the representations of 
SU(3)c x SU(3)l x SU(3)r as follows: 

L~(l,3,3), Ql ~ (3, 3,1), Qr ~ (3,1,3). (2) 

In matrix notation, L is a 3 x 3 matrix, Qr is a column 
vector, and Qr is a row vector: 


( L\ 

E~ 

e \ 


( u b \ 


L = ( E+ 

L | 


, Q b L 

= d b , 


\ e+ 

V 

Ll) 


\D b J 





Qr 

= (u b d b D b ). 

(3) 

b = 1,2,3 

is a 

color index. 

The components 

u, d 


are the left-handed up-type and down-type quarks from 
the SM, whereas I? is a new quark with electromagnetic 
charge —The components u, d, D are their respective 
right-handed counterparts. The lepton field contains the 
charged leptons and the left-handed neutrino v. It 
contains several new states: a right-handed neutrino D; 
three neutral states L\, L\, L 3 ; and a pair of charged lep¬ 
tons E±. A generation index a = 1,2, 3 on the fermion 
fields in eq. (3) has been suppressed. 

The Higgs fields Hi, H 2 cannot both couple to 
fermions, since this would lead to flavor-changing neutral 
current (FCNC) processes, which are severly restricted 
by experiment. In order to suppress tree-level FCNC in¬ 
teractions, the existence of a Z 2 -symmetry is assumed 
under which Hi (H 2 ) is even (odd). The fermions are 
even under this symmetry as well, which implies that H 2 
does not couple to fermions. The Yukawa couplings are 
of the form 

Cy = - gtG aP (qZHTqZ + \f k e lmn L\V m (Hi) k ^j 

- A aP (Q a R H T Aq Q p L + e^ k L\U m (H Al )\ lm} ) 

- vb (G2)a/3Tr {^} Tr W L ' 3 } + h - c - 

(4) 


The first line is a Yukawa interaction built from the fields 
we have already introduced: the parameter g t is a di¬ 
mensionless coupling, G a p is a symmetric 3x3 gener¬ 
ation matrix, and e is a totally antisymmetric symbol 
with e 123 = £123 = +1. This interaction is sufficient 
to reproduce the up-quark masses by choosing a gen¬ 
eration basis in which gtG a p is diagonal and fitting its 
diagonal components to the up-quark masses [23]. The 
vev Mi in eq. (1) gives large masses to D, E ± . The 
second line in eq. (4) contains interactions with new 
scalar fields H Aq ~ (1,3,3) and Hai ~ (1,3,6), cou¬ 
pling to the fermions with a Hermitian antisymmetric 
matrix A a f 3 .[28] These terms come from the couplings of 
a scalar field H A in the antisymmetric 351a represen¬ 
tation of Eq. This interaction is necessary to describe 
the masses and mixings of the down quarks and charged 
leptons correctly: the Standard-Model down quarks and 
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charged leptons are mixed with their heavy partners via 
the seesaw mechanism. A good fit for the masses and 
mixings of the Standard-Model fermions is obtained us¬ 
ing only very few extra parameters [23]. It is assumed 
that the fields Ha q , Hai have negligible mixing with Hi, 
H 2 in order to simplify the analysis of the scalar spec¬ 
trum. 

At this stage, neutrinos are still Dirac particles with 
masses comparable to the other fermion masses. The 
third line in eq. (4) is necessary to obtain neutrino 
masses in accordance with experiment. This line con¬ 
tains an effective dimension-five Yukawa interaction that 
could originate from the exchange of a new heavy Dirac 
fermion that is a trinification singlet [20]. It violates 
the ^-symmetry and mixes the neutrinos v, v with the 
other neutral leptons L \, L 2 , L\, giving rise to a gener¬ 
alised seesaw mechanism. The light-neutrino mass ma¬ 
trix introduces two additional parameters, which can be 
fixed by the experimentally observed atmospheric mass- 
squared difference and the lightest neutrino mass [23]. 


III. AN EFFECTIVE TRINIFICATION MODEL: 
THE LET MODEL 


We now consider the effective trinification model ob¬ 
tained after integrating out the Higgs fields that are made 
heavy by the large vevs M\ and M 2 . We expect the Higgs 
fields that are right-handed singlets with respect to the 
SM to obtain large masses. These are the fields with the 
SU(3)r index (3). On the other hand, the fields (3,1) and 
(3,2) of H 2 that are left-handed singlets (with respect to 
the SM) are kept. These Higgs fields are necessary to 
describe the breaking of the left-right symmetry. The 
corresponding scale M is certainly larger than the weak 
scale, but may be much lower than Mi and M 2 - 

The 2x2 blocks in the upper left corners of the fields 
Hi, H 2 transform as bidoublets $i ,$ 2 ~ (1,2, 2,0) un¬ 
der SU(3) C x SU{2) l x SU{2) r x U{1) B -l- The (3, 1 ) 
and (3, 2) components of H 2 transform as a right-handed 
doublet ~ ( 1 , 1 , 2 , 1 ): 




(®i, 11 ®i,2l\ 

U+12 




** = (*£ K) ** 


(m\ ml o\ 
mi mi 0 , 

\ 0 0 0 / 

/ 0 0 0 \ 

0 0 0 

\m)l { h 2 )1 0 ] 


(5) 


These fields obtain the following vevs: 




1 fvj 0 \ 


(®r) 


> M >- 


( 6 ) 


The fields $ 1 , $ 2 , &r, and their vevs are sufficient to de¬ 
scribe the symmetry breaking from SU(3)c x SU(2) L x 
SU (2 )rX U (1 )b-l to electromagnetism via the SM. Note 
that a right-handed doublet like resides in Hi as well. 


In principle, $ 1 ? may be a combination of both. Note that 
the (3,3) components of Hi and H 2 are total gauge sin¬ 
glets. A combination of them, if sufficiently light, could 
be a dark matter candidate. Also note that $2 and <f>u 
do not couple to fermions. 

Besides the eight gluon fields of SU(3)c, the gauge- 
boson sector consists of seven vector gauge bosons. Their 
masses and mixings are given in appendix A. The vevs 
of $ 1,2 contribute to the W mass and are therefore 
constrained by the relation vf + b\ + v 2 + b\ = v 2 = 
(246 GeV) 2 . It is convenient to reparanreterize the vev 
parameters as 

Vi =v cos a cos (3i, V 2 = fsinacos/ 32 , 

bi =i> cos cc sin/?i, 62 = v sin a sin /3 2 ■ (7) 

The parameter M is a mass scale above the electroweak 
scale. Because of this hierarchy, it will often be conve¬ 
nient to use the small parameter £ = v/M. 

We build the scalar potential of the LET model from all 
possible gauge-invariant renormalizable operators con¬ 
sisting of $ 1 , <b 2 , leaving out those that could not 
have arisen from the trinification model. This means 
that we do not include operators involving charge con¬ 
jugates: since the 3 and 3 representations of SU (3) are 
inequivalent, such operators have no possible origin in 
the trinification model. The resulting potential is 


V =Hi($!, $ fl ) + H 2 ($i, $ 2 ,*«), 

VI =yTr {$!$!${$!} + ^(<b J? .^) 2 

+ A 4 Tr{$j<bi} ($R$k) + A 5 $R<bI$i$k 
+ huTr {$!$1 j + Ar®R®r + (hi det 4>i + h.c. ) , 
H 2 =yTr {$|$ 2 } 2 + yTr{$t$ 2 $|$ 2 } 


+ A 3 T 1 - 


{$2^2} ($«$],) + A 4 <f> fi <f 4 <f> 2 $R 
A 5 TV<j$l$i}Tr{$!$ 2 } +A 6 iTrjV^}' 

^ {<I>j'I> 2 } 

A 


Tr- 


+ h.c. ) +A 8 Tr{$l$i$^ 2 } 


+ A 9 Tr + ~y (Tr + h.c. 

+ h 22 Tr {$ 2 ^ 2 } + (hi det $ 2 + h.c. ) . ( 8 ) 


For later convenience, we have split the potential into 
the parts Vj and V 2 that respectively do and do not de¬ 
pend on <f> 2 . As mentioned in the introduction, the five 
dimensionful parameters hu, h 22 : hij: hi 1 A are n0 free 
parameters. They are determined in terms of the vev pa¬ 
rameters Vi, bi, V 2 , b 2 , M and the dimensionless param¬ 
eters A i, A j by the requirement that V has an extremum 
at the appropriate place. 

In the trinification model, the fermions obtain their 
masses from Yukawa interactions with H 4 . In order to 
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combine the leptons and $1 into a gauge singlet, we need 
the antisymmetric tensor ia^. which can be absorbed into 
a redefinition of the lepton fields. Absorbing a minus sign 
into the phase of the e ± fields, the fermionic field content 
of the LET model becomes 

Ql = (fj ~ (3,2,1, |), Qr = (u d) ~ (3,1,2,-i), 

L~= ^ -(1,2,1,-1), L+ = (u e+) -(1,1, 2,1). 

(9) 

After integrating out the heavy fields, the first line of the 
Yukawa Lagrangian in eq. (4) becomes 

C Y = -G a p (Q a R ^>lQ p l + L +a <f>rL~^ + h.c. (10) 

This term is insufficient to describe all fermion masses 
correctly: down-quark masses as well as the charged- 
lepton and neutrino masses would be proportional to each 
other and the CKM matrix would be a unit matrix at this 
point. In order to describe the fermion masses and mix¬ 
ings correctly, interactions with the additional fields ip 4 g , 
Hai need to be included. Here we restrict ourselves to the 
single Yukawa term in eq. (10) and fit the free parame¬ 
ters to the top- and bottonr-quark masses: these fermions 
are the most relevant to compare our analysis to experi¬ 
mental searches for new physics. We assume that flavor 
physics for the lower fermion masses does not influence 
the spectrum of the scalar particles. Fitting the free pa¬ 
rameters to the top- and bottom-quark masses, we find 
v = 246 GeV, tan/3i = m,b/in t => f3\ = 0.0166. Since <f> 2 
does not contribute to the fermion masses, we have no 
such restrictions on /3 2 and a. 


A. The Single-Bidoublet model 

The scalar potential of the LET model in eq. (8) con¬ 
tains 15 free coupling constants, which are difficult to 
deal with. As an intermediate step towards an under¬ 
standing of the LET model, we first discuss a model 
from which <I> 2 has been omitted; we will refer to this 
setup as the Single-Bidoublet (SB) model. This model 
corresponds to the LET-model limit a = 0, ^2 oo 
(a = 0 implies V 2 = &2 = 0, for which /r| 2 is n o longer 
constrained by the location of the minimum of the scalar 
potential). Note however that the SB model is not an ap¬ 
propriate effective field theory of the trinification model, 
since the vev parameters of <I> 2 are of order v for the gen¬ 
eral case a ^ 0. Rather, we use it as a toy model that 
helps us study the phenomenology of the LET model. 

The most general scalar potential for the SB model is 
given by Vi(<f>i, <b/j) in eq. (8). The minimalisation of 
this potential at the vev in eq. (6) fixes the dimensionful 
parameters /t 2 -^ /t^, /r 2 in terms of the \t and the vev 
parameters V \, b\, M. The vev parameters iq, bi can be 
reparameterized in terms of v, /3i by eq. (7), with a = 0. 


The scalar fields <f>i, &r contain twelve real scalar com¬ 
ponents in total. After spontaneous symmetry breaking, 
six of them become massless Goldstone bosons that give 
mass to the six massive vector gauge bosons. The remain¬ 
ing components mix to form six massive scalars: three 
CP-even scalars /i°, P°, and H®, one CP-odd scalar A 0 , 
and a pair of charged scalars H ± . Their definitions in 
terms of gauge eigenstates as well as their masses are 
given in appendix B. We identify hr with the Standard- 
Model-like Higgs particle that has been observed at the 
LHC [1, 2], since it is the only scalar that naturally has 
a mass at the electroweak scale. The other scalars have 
masses of order M unless we give O (£ 2 ) = O ( v 2 /M 2 ) 
values to some of the dimensionless scalar-potential pa¬ 
rameters. 

Note that the SB model resembles the two-Higgs- 
doublet model (2HDM) [27] (see [29] for a recent review). 
The 2HDM is an extension of the SM in which the scalar 
sector contains an additional SU(2)l doublet. It has 
been studied extensively since it provides a low-energy 
description of various models such as supersymmetry (see 
e.g. [30] for a review), composite Higgs models [31], and 
little Higgs models [32]. 

It is easy to see why the SB model resembles the 2HDM 
at the Lagrangian level. To this end, we write 4>i = 
(ic^i,<(> 2 ) and = (0+, 0o), where 0 1, 2 are SU(2) L 
doublets and 0 +i 0 are SU{2)l singlets. Since $r has a 
much larger vev M than the vev components tq, b\ of 
'Ll, the mixing among $ 1 , will be of order 1. If 
we set = 0 in the scalar potential V\ in eq. (8), we 
can rewrite the entire scalar potential in terms of 0i, 0 2 
only. The result is a 2HDM potential (see e.g. eq. (98) 
in ref. [29]) with the following constraints on the scalar 
parameters: 


\ 2HDM 
a 1,2,3 


—Ai + A 2 , 


m ll — m 22 —Mil, 


\2HDM _ \ 

+4 — ~+2, 

2 

TO 12 = -H V 


\2HDM _ n 
A 5,6,7 — U ’ 


( 11 ) 


We can rewrite the Yukawa sector of the LET model in 
terms of 0i, 0 2 as well. In the Lagrangian in eq. (10), the 
up-type fermions couple only to 0 1 whereas the down- 
type fermions couple only to 0 2 . Hence the SB model 
resembles a constrained type-II 2HDM setup. It differs 
from the 2HDM due to the presence of two additional 
heavy SU{2)r singlets. The neutral one gives rise to an 
additional physical particle P® !h a ! is fermiophobic. Its 
mass can be tuned independently from the masses of the 
2HDM-like scalars P°, A °, (see eq. (B2)). This can 

result in phenomenological scenarios that cannot appear 
in the 2HDM. Moreover, the vev ratio tan/3 = (02)/(0i) 
is a free parameter in the 2HDM, whereas tan/3i = 
rrib/mt is fixed in the LET model. Likewise, the mix¬ 
ing angle a 2HDM of the CP-even scalars in the 2HDM 
can be taken to be a free parameter: the A 2HDM are usu¬ 
ally rewritten in terms of the scalar masses and a 2HDM . 
However, the three mixing angles of the CP-even scalars 
in the SB model cannot be treated as free parameters: 
they are approximately fixed by the value of /3i unless at 
least one of the new scalars becomes light. 
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B. LET-model scalar spectrum 

Now let us consider the scalar sector of the LET model. 
Compared to the SB model, it contains an additional 
bidoublet $2 with eight real scalar components. This 
makes 14 physical scalars in total: five CP-even states 
h°. Hi, HH®, H%; three CP-odcl states AA®, 
and three pairs of charged states H±, H±, . Their 

definitions and masses are given in appendix B. The CP- 
even states h°, H® have masses of order v, whereas H 
H®, H 4 have masses of order M. This is not surprising: 
if we decouple the bidoublet $ 2 from the model, we get 
one light state h° and two heavy states. Since $1 and <f> 2 
are copies of the same representation, we expect that $ 2 
adds one light and one heavy scalar to the spectrum as 
well. The CP-odd state Ai is light, whereas A°, ^3 are 
heavy. Again, this is not surprising. In the SB model, 
the CP-odd components of the bidoublet $1 give rise 
to one Goldstone and one heavy state. Thus we would 
expect $ 2 to contribute one heavy state as well. Since 
there are no more would-be Goldstones, the other CP- 
odd component of <b 2 becomes a massive state with a 
mass of order v. Similarly, is light whereas H±, H± 
have masses of order M. 


IV. TRINIFICATION PHENOMENOLOGY 

A. The SB model 

Now we turn to the phenomenological scenarios that 
are allowed by the scalar sector of the trinification model. 
As a first step towards understanding the trinification 
phenomenology, we consider the phenomenological sce¬ 
narios that are possible in the SB model. The free pa¬ 
rameter space of the SB model is spanned by M and the 
five scalar parameters A*. A full analysis of this parame¬ 
ter space and the possible signatures is beyond the scope 
of this work. Instead, we define a set of benchmark points 
that lead to distinct phenomenological features. To get a 
feel for the possibilities, consider the scalar masses given 
in eq. (B2). The mass of h° can be adjusted by changing 
the values of Ai, A 2 , A3, A4. We tune these parame¬ 
ters such that m/jO = 126 GeV for each benchmark. The 
leading contributions to m H 0, m^o, m#± are all given by 
yj A5 M , so we expect them to have similar masses, with 
O (v) mass splittings. On the other hand, m H a is pro¬ 
portional to y/X^M , which can be tuned independently of 
the other scalar masses. Thus we expect the SB model to 
allow for different mass hierarchies or compressed spec¬ 
tra, depending on the magnitudes of A3, A5. 

If any of the parameters A3, A5 have O (£ 2 ) values, 
some of the new scalars may obtain O (y) masses. Thus 
a number of phenomenologically different scenarios are 
possible. If A 3 , A 5 are not too small, all new scalars 
obtain O ( M ) masses, beyond experimental reach. We 
denote these benchmarks as having a Single Large Hier¬ 


archy (SLH); an additional hierarchy between H /, A 0 , 
H± on one side and on the other side is possible, but 
of little phenomenological interest. For A 5 ~ O (£ 2 ) but 
not too small A 3 , the 2HDM-like scalars H®, A 0 , H ± all 
have O (i>) masses and could be observed at the LHC. We 
choose two parameter sets such that these particles have 
masses in the O (100 GeV) ballpark, and denote these 
benchmarks as 2HDM-1 and 2HDM-2 since they have a 
2HDM-like spectrum at low energies. For A 3 ~ O (£ 2 ) 
but sizable A 5 , the fermiophobic state H® lies within ex¬ 
perimental reach[33]. We choose two parameter sets such 
that m H 0 lies in the O (100 GeV) range and denote these 
benchmarks as Light Fermiophobic (LF). Note that for 
small A 3 , the parameters A 4 , A 5 need to be chosen suf¬ 
ficiently small to ensure that ?nj) 0 remains positive (see 
eq. (B2)). This means that the 2HDM-like scalars have 
masses well below M. A combination of the 2HDM-like 
and LF scenarios is possible as well: if both A 3 and A 5 
are sufficiently small, all new scalars could have masses 
within experimental reach. Again, we choose two param¬ 
eter sets and refer to these benchmarks as Compressed for 
having a spectrum compressed around the electroweak 
scale. 

We also consider two special cases of the LF sce¬ 
nario. For small enough A 3 , TL® can h av ® a mass in the 
O (1 GeV) ballpark. Such a state would decay into pairs 
of photons only, since H% is fermiophobic. If the sig¬ 
nal strength for its decay is low enough, it could have 
escaped detection so far; we discuss the relevant exper¬ 
imental constraints in more detail in section VI. We re¬ 
fer to this scenario as having a Very Light Fermiophobic 
(VLF) Higgs. The second special case is when h° and H® 
are approximately degenerate. If their mass difference is 
less than their widths, both states would contribute to 
the signal strength used for the Higgs discovery, lead¬ 
ing to a ‘twin Higgs’[34] [35, 36]. This might result in 
deviations of the measured Higgs couplings from their 
Standard-Model values. We tweak the parameters such 
that m h o jjo = 126 GeV and denote the corresponding 
benchmarks as Twin-1 and Twin-2. 

For each benchmark point, we use the values v = 246 
GeV, /?i = 0.0166 as well as the experimental values 
sin 2 9w = 0.23126, c/l = 0.65170, and our best fit 
9' w = 0.62 as in appendix A. The other gauge cou¬ 
plings are fixed by the identities gn = gL tan 0w/ sin 9' w , 
2 g' = g l tan 9\\r j cosf V w . As for the scale M, we consider 
both a high scale M = 10 10 GeV well outside experi¬ 
mental reach and a lower scale M = 10 4 GeV just be¬ 
yond LHC reach. We ensure that the constraints for vac¬ 
uum stability and S-matrix unitarity are satisfied. Multi- 
Higgs potentials in general can have several minima. De¬ 
termining the global minimum of the potential is already 
challenging in the 2HDM [37]. We could not check that 
our minimum is the global one. The corresponding pa¬ 
rameter values are given in table I. Note that except for 
the SLH scenario, all benchmarks contain very small val¬ 
ues for A 3 , A 4 and/or A 5 : these parameters need to have 
O ( v 2 /M 2 ) values to compensate for the large vev of the 
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Benchmark point 

m H o 

m^o 

m H ± 


SLH-1 

3.2 • 10 9 

3.2 • 10 9 

3.2 • 10 9 

6.9 • 10 9 

SLH-2 

3.2 • 10 3 

3.2 • 10 3 

3.2 • 10 3 

6.9 • 10 3 

2HDM-1 

488 

488 

m 

6.6 ■ 10 9 

2HDM-2 

488 

488 

m 

6.6 • 10 3 

LF-1 

3.9 • 10 6 

3.9 • 10 6 

3.9 • 10 6 


LF-2 

5.5 • 10 3 

5.5 • 10 3 

5.5 • 10 3 


Compressed-1 

496 

496 



Compressed-2 

496 

496 



VLF-1 

1.0 • 10 6 

1.0 • 10 6 

1.0 • 10 6 


VLF-2 

5.9 • 10 3 

5.9 • 10 3 

5.9 • 10 3 


Twin-1 

7.1 • 10 4 

7.1 • 10 4 

7.1 • 10 4 


Twin-2 

2.2 • 10 3 

2.2 • 10 3 

2.2 • 10 3 



Benchmark 

M 

Ai 

^2 

^3 

a 4 

A 5 

SLH-1 

10 10 

0.24 

0.24 

0.47 

0.32 

0.2 

SLH-2 

10 4 

0.24 

0.24 

0.47 

0.32 

0.2 

2HDM-1 

10 10 

0.41 

0.4 

0.44 

0.49 

5 • 10" 15 

2HDM-2 

10 4 

0.41 

0.4 

0.44 

0.49 

5•10“ 3 

LF-1 

10 10 

0.133 

0.13 

2 • 10~ 15 

1 • 10" 12 

3•10“ 7 

LF-2 

10 4 

0.14 

0.14 

2•10 -3 

5.5 ■ 10" 3 

0.6 

Compressed-1 

10 10 

0.133 

0.13 

1.1 • 10” 15 

1 • 10" 12 

5 • 10 -15 

Compressed-2 

10 4 

0.15 

0.14 

1.1 ■ 10 -3 

5.2 ■ 10 -3 

5-10“ 3 

VLF-1 

10 10 

0.133 

0.13 

1 • 10~ 2 ° 

1 • 10" 14 

2•10“ 8 

VLF-2 

10 4 

0.13 

0.13 

4.5 • 10“ 7 

4•10- 5 

0.7 

Twin-1 

10 10 

0.13 

0.133 

1.58 • 10 -16 

1 • 10" 16 

1 • 10" 10 

Twin-2 

10 4 

0.131 

0.131 

1.59 • 10~ 4 

1•10^ 5 

0.1 


TABLE I. Definitions of the benchmark points for the SB 
model in terms of the free parameters M (in GeV) and A;. 
The parameter values v = 246 GeV, /3i = 0.0166, sin 2 9w = 
0.23126, gL = 0.65170, 9' w = 0.62, g R = gLt&ndw/ sia9' w , 
2 g' = gL tan 9w / cos 9\y are kept fixed. 

corresponding scalar invariants. This makes these sce¬ 
narios unnatural. 

To calculate the scalar masses, we do not use the ap¬ 
proximations in eq. (B2), since subleading terms in £ may 
become large for small A 3 and/or A 5 . Instead, we eval¬ 
uate the mass matrix numerically in Mathematica [38] 
and then extract the masses and mass eigenstates. The 
corresponding particle masses are given in table II. The 
scalar mass eigenstates are almost equal to the gauge 
eigenstates in most benchmarks: mixings are mostly be¬ 
low the percent level. In the Compressed-2 benchmark, 
there is a 2% mixing of h\ n with h^. Only Twin-2 gives 
large scalar mixing: in terms of squares of amplitudes, 
the SM-like Higgs h° is 62% n and 38% h ° R , whereas 
the fermiophobic Higgs Lf® is 38% h\ n and 62% h R ; 
mixing with h\ 22 is negligible. 

B. The LET model 

In the SB model, the new physics decouples from the 
SM unless some of the dimensionless scalar parameters 
are set to O (£ 2 ) values. The reason is the fact that there 
is a large hierarchy v <C M among the vev parameters 
of $!, <f>^. The picture changes with the inclusion of the 
bidoublet $ 2 , since its vev parameters are bounded by 
the electroweak scale. In the absence of a large hierar¬ 
chy between the vev components of $1 and $ 2 , we ex¬ 
pect that the LET model allows for large mixing between 
the components of both bidoublets if the dimensionless 
scalar-potential parameters have O (1) values. Hence the 
model naturally contains new scalar particles with O (u) 
masses. As in the SB model, significant mixing of $ 1 , $2 
with the components of <!># is only expected for unnatu¬ 
rally small values of some of the parameters. 


TABLE II. Scalar masses for the benchmark points defined in 
table I. All masses are given in GeV. The mass of the SM-like 
Higgs h° has been tuned to 126 GeV in each case. 


These insights have important consequences for the 
benchmark scenarios considered in the previous section. 
The SLH scenarios have no analogon in the LET model: 
there is always new physics within experimental reach. 
The 2HDM-like, LF, VLF and Twin scenarios become 
natural possibilities, since the 2HDM-like scalars can nat¬ 
urally have O ( v ) masses. Moreover, mixing of $1 and $2 
may make these scalars fermiophobic, in which case this 
scenario could be distinguished from the usual 2HDM. 
Hence we expect the LET model to be predictive: it al¬ 
lows for phenomenologically interesting, experimentally 
testable scenarios, since the new physics does not decou¬ 
ple from the SM in the large-M limit. 

In the following, we discuss the prospects for measur¬ 
ing the Higgs-coupling modifications in the context of the 
LET model. To this end, we define a new set of bench¬ 
mark points, inspired by the considerations given above. 
As a starting point for choosing the parameter values, we 
observe the following about the scalar masses in eq. (B3). 
The main contribution to m/jO is given by Ai + A 2 cj^, with 
an overall scaling factor (? a due to the presence of the sec¬ 
ond bidoublet. This means that a smaller a is generally 
accompanied by a smaller value for Ai + . Similarly, 

m H 0 is mainly determined by Ai + A 2 c ^ 2 with an overall 
factor s 2 a , so smaller values of a should be compensated 
by larger values of Ai+A 2 C ^ 2 . The h°—H° mass difference 
and mixing are governed by S( 2a ) as well as the scalar pa¬ 
rameters A5 i 6,7,8,9,io> so we tune these parameters until 
we have a parameter set that corresponds to the desired 
benchmark scenario. The scalar H ° has m H 0 ss 

so we should take A 3 > 0. The squared masses of A°, 
are determined by A 6 , 7 , 9 ,io with an overall minus sign, so 
we take these parameters to be negative to guarantee a 
positive-definite mass matrix. Positivity of the squared 
masses of H 2 , A 2 , H 2 requires A 5 > 0, and the squared 
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2HDM-3 

2HDM-4 

VLF-3 

Twin-3 

Twin-4 

sin a 

0.93 

0.43 

0.50 

0.72 

0.33 

sin /?2 

0.17 

0.65 

0.17 

0.16 

0.11 

Ai 

1.0 

0.20 

0.13 

0.34 

0.17 

A 2 

1.0 

0.18 

0.13 

0.35 

0.16 

^3 

0.50 

0.50 

0.50 

0.50 

0.42 

a 4 

0.010 

0.12 

0.13 

0.27 

0.12 

A5 

0.20 

0.50 

0.20 

0.20 

0.50 

Ai 

0.40 

1.3 

0.27 

0.34 

1.3 

A2 

0.40 

1.3 

0.27 

0.34 

1.2 

A3 

0.27 

0.20 

0.27 

0.27 

0.19 

a 4 

0.20 

0.10 

0.20 

0.20 

1.0 

A5 

0.20 

0.046 

0.54 

0.32 

0.060 

A 6 

-0.40 

-0.50 

-0.43 

-0.42 

-0.30 

A7 

-0.24 

-0.20 

-0.24 

-0.24 

-0.30 

As 

0.84 

0.95 

0.83 

0.84 

1.0 

A9 

-0.050 

-0.10 

-0.04 

-0.053 

-0.10 

A10 

-0.30 

-0.30 

-0.30 

-0.30 

-0.30 


TABLE III. Benchmark-point definitions for the LET model. 
The values v = 246 GeV, M = 10 10 GeV, di = 0.0166 are 
kept fixed. 



2HDM-3 

2HDM-4 

VLF-3 

Twin-3 

Twin-4 

m h o 

TTlrrO 
n 1 

rrirro 
n 2 

TTlrrO 
n 3 

TTlrrO 

126 

182 

3.2 ■ 10 9 

3.3 • 10 9 
7.1 ■ 10 9 

126 

148 

5.0 • 10 9 
5.8 • 10 9 
7.1 • 10 9 

126 

3.9 

3.2 • 10 9 

3.3 • 10 9 
7.1 • 10 9 

126 

126 

3.2 • 10 9 
3.2 • 10 9 
7.1 • 10 9 

126 

126 

5.0 • 10 9 
7.2 • 10 9 

6.5 • 10 9 

m A 0 

m A 0 

m A 0 

179 

3.2 ■ 10 9 

3.3 • 10 9 

134 

5.0 ■ 10 9 
5.8 • 10 9 

179 

3.2 • 10 9 

3.3 • 10 9 

179 

3.2 • 10 9 
3.2 • 10 9 

190 

5.0 • 10 9 
7.2 • 10 9 

m H ± 

TTlrr± 

H 2 

TTlrr± 

__ 

171 

3.2 ■ 10 9 

3.3 • 10 9 

135 

5.0 • 10 9 

5.8 • 10 9 

173 

3.2 • 10 9 

3.3 • 10 9 

173 

3.2 • 10 9 

3.2 • 10 9 

173 

5.0 • 10 9 

7.2 • 10 9 


TABLE IV. Scalar masses for each of the benchmark points 
defined in table III. All masses are in GeV. 

masses of H®, Ag, require A 4 /c( 2 / 3 2 ) >0. We tune the 
parameters such that = 126 GeV is fixed and enforce 
a set of vacuum stability conditions. The parameter val¬ 
ues for the benchmark points are given in table III, and 
the resulting scalar masses are given in table IV. 

We have defined two 2HDM-like benchmarks, in 
which the scalars H Ag, H ^ all have masses in the 
O (100 GeV) range. In 2HDM-3, the state h° is almost 
purely h° n , whereas H® is 97% h 211 and 3% h 2 22 - The 
state A° (H±) is 87% a® n (hf 2 1 ) and 13% a 2 n {hf 2 1 ). 


The mixings are even larger in the 2HDM-4 benchmark: 
h° is 74% hi tll , 15% h 211 , and 11% h 2 22 whereas H° is 
25% /i? n , 43% /i 2 n and 32% h 2 22 . The state Ag (H±) 
is 18% a,i n {hf 21 ), 47% a 211 (h 221 ), and 35% a 222 
(/i^i 2 ). In both cases, the 2HDM-like scalars have sup¬ 
pressed couplings to fermions with respect to the type-II 
2HDM. 

We also define the VLF-3 scenario, in which H® has 
a mass of only a few GeV. The states Ag, have 
masses within experimental reach. The scalar mixing 
is significant: h° is 65% 34% h 211 , and 1% h 2 22 , 

so we expect reduced fermion couplings. The state H® 
is 35% /i5,u> 63% ft®,ii) and 2% h 2 22 . Since it is light 
and mainly fermiophobic, it could have evaded the LEP 
searches. The state Ag is 25% Og n (h^ 21 ), 73% 

a 2 ,n (^ 2 , 21 )' and 2% a 2 22 (^ 2 , 12 )' 

Furthermore, we define two Twin benchmarks with dif¬ 
ferent amounts of scalar mixing. For Twin-3, the state 
h° is 87% hi U and 13% h 2 ll , whereas H° is 13% h^gg, 
85% h 211 , and 2% h 2 22 . The lightest CP-odd and 
charged states are almost 50-50 mixtures of fermiophilic 
and fermiophobic states: A° (H±) is 51% a? n (hf 21 ), 
47% a 211 (hf 21 ), and 1% a 222 (^ 2 ) 12 )- For Twin-4, h° 
is almost purely h\ 13i whereas H ° is 99% h 211 and 1% 
h 2 22 ■ The lightest CP-odd and charged states are mostly 
fermiophobic: A? (H±) is 11% a° )11 {hf 21 ), 88% a 211 
(^ 2 , 2 l)j anc l 1% a 2,22 (^ 2 , 12 )’ 


V. HIGGS-COUPLING MODIFICATIONS 

In the SM, the Higgs couplings are fixed in terms of the 
particle masses and the vev of the Higgs field. Hence an 
independent measurement of these couplings provides an 
important test of the SM. These couplings are generally 
modified in the presence of an extended Higgs sector [39]. 
The Higgs-coupling modifications A x are defined as the 
deviations of the Higgs couplings g x = gh° X x from their 
SM values, where x is any SM particle and gh° xx is the 
coefficient of the operator h°xx in the Lagrangian: 

g x = (1 + A x )g* M . (12) 

The loop-induced Higgs coupling to photons can be writ¬ 
ten as follows: 

57 = (1 + A 7 M + A 7)57 M ( 13 ) 

Here A® M is the coupling modification that is induced by 
coupling modifications of the Standard-Model particles 
generating the coupling. The term A 7 represents contri¬ 
butions from non-SM particles running in the loops. The 
Higgs-coupling modifications have been extracted from 
LHC data using the tool SFitter [39-42] (see fig. 1). 






FIG. 1. Most recent fit of the Higgs-coupling modifications 
to LHC data. The red points correspond to the expected 
SM result A x = 0, whereas the dark blue points give the 
results from the data if the photon coupling is assumed to 
be determined by the W and t loops only. The light blue 
points give the results if a free coupling shift in A 7 due to 
new physics is allowed. Figure taken from ref. [39]. 


A. The SB model 

The Higgs-coupling modifications of the SB model are 
given in appendix C. In order to see whether our bench¬ 
mark scenarios could be distinguished from the SM in ex¬ 
periment, we calculate the Higgs-coupling modifications 
numerically for each benchmark point. The results are 
listed in table V. 

For the SLH benchmarks, the modification of the quar- 
tic Higgs self-coupling is as large as 83%, whereas the 
other couplings have negligible deviations from their SM 
values. Hence if we could measure both Higgs self¬ 
couplings, a SB model with a large hierarchy could be 
distinguished from the SM by the strength of the quar- 
tic Higgs self-coupling, even if M is large. Like in the 
2HDM [43], the self-couplings are generally modified due 
to the more complicated structure of the scalar potential, 
although these modifications may also vanish along some 
directions of the parameter space (see appendix C). 

For the 2HDM-like benchmarks, there is an 11% in¬ 
crease of the Higgs coupling to b quarks. The reason 
is that the main contribution to Ain this scenario is 
proportional to £ 2 /A 5 ~ 0.1. The W, Z and t couplings 
obtain similar modifications, but these are suppressed by 
factors ~ 10~ 3 and ~ 10 -4 respectively. Hence 
the 2HDM-like hierarchy is characterized by an increase 
in only the b coupling. The quartic Higgs self-coupling is 
enhanced by a factor 3. 

The LF, Compressed, and VLF benchmarks have cou¬ 
pling modifications at or below the percent level. These 
scenarios would not be distinguishable from the SM via 
the Higgs-coupling modifications. The Twin-1 bench¬ 
mark has very small coupling modifications as well, but 


the Twin-2 benchmark shows 20% reductions of all tree- 
level Higgs couplings and 50-60% decreases of both self¬ 
couplings. Currently, the errors on the Higgs couplings 
are still large enough to allow a 20% deviation (see fig. 1). 
If the errors can be reduced after the 14 TeV run, the 
Twin-2 benchmark can be put to the test. 


B. The LET model 

Because of the large number of free parameters and 
scalar mixing angles, we restrict ourselves to a numeri¬ 
cal analysis of the Higgs-coupling modifications for the 
benchmark points defined in section IV B. The result¬ 
ing coupling modifications have been summarized in ta¬ 
ble VI. Note that contrary to the SB model, the photon¬ 
coupling modification is not negligibly small any more: 
the charged scalar has a mass of order v and hence 
yields a sizable contribution to the effective photon cou¬ 
pling. 

The 2HDM-3 benchmark has large coupling modifica¬ 
tions: the couplings to V = W,Z are suppressed by a 
factor 0.3, whereas the couplings to t, b are enhanced 
by almost a factor 3. This is not surprising: the W, Z 
couplings of the SB model are proportional to v, whereas 
the t, b couplings are proportional to m^b/v. In the LET 
model, we have to substitute v vcosa (see eq. (7)). 
Since h° is almost purely h\ n in the 2HDM-3 bench¬ 
mark, the W, Z couplings are suppressed by cos a = 0.36 
whereas the t, b couplings are enhanced by 1/cos a = 2.8. 
This benchmark point is clearly incompatible with the 
measured coupling modifications in fig. 1. In contrast, 
the 2HDM-4 benchmark has smaller but still sizable cou¬ 
pling modifications. The V couplings are suppressed by 
cos a = 0.90, but the total coupling is a few percent 
higher because of the contribution from the second bidou¬ 
blet. On the other hand, the quark couplings are en¬ 
hanced by a factor 1/ cos a = 1.11, but the total coupling 
modifications are negative: since h° contains a significant 
admixture of the fermiophobic $ 2 , the t. b couplings are 
reduced. The coupling modifications for 2HDM-4 are 
consistent with the measured coupling modifications. 

The VLF-3 scenario has large mixing between the 
fermiophilic and fermiophobic scalar gauge eigenstates. 
As such, there is again a tension between cos a and 
$1 — 4>2 mixing. The resulting coupling modifications 
are at the percent level, all compatible with the mea¬ 
sured values. 

Like the 2HDM-3 benchmark, the Twin-3 scenario has 
large coupling modifications. The W, Z couplings are 
reduced by about 60%, mostly due to interference be¬ 
tween the contributions of the fermiophobic and fermio¬ 
philic scalar components. The t, b couplings are enhanced 
by 34%: the factor 1/ cos a = 1.4 is slightly reduced by 
$1 — 4>2 mixing. This benchmark point is incompati¬ 
ble with the data. However, the Twin-4 scenario has 
percent-level coupling modifications, compatible with the 
measured values. This is because the V couplings are 
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Benchmark point 

A 

w 

A z 

A 

t 

A b 

A 7 

A a 

3 h 

A *4h 

SLH-1 

-1.4 • 

10— 16 

0 

0 

-1.4 • 

10~ 16 

0.0 

-5.0 • 10“ 17 

0. 

0 

0.83 

SLH-2 

-1.4 • 

10^ 4 

-4.1 

10“ 4 

-1.4 • 

10~ 4 

2.1 • 10~ 3 

-5.0- 10“ 5 

-4.2 • 

10~ 4 

0.83 

2HDM-1 

-1.6- 

10^ 6 

-1.6 

10“ 6 

-3.2 • 

10~ 5 

0.11 

-1.7- 10 -3 

-9.7- 

10~ 5 

2.1 

2HDM-2 

-3.8- 

10- 4 

-6.4 

10“ 4 

-4.1 • 

10~ 4 

0.11 

-1.6- 10“ 3 

-1.2 • 

10~ 3 

2.1 

LF-1 

-6.2 • 

10^ 7 

-6.2 

10- 7 

-6.2 • 

10- 7 

-6.2 • 10' 7 

-2.6 • 10“ n 

-1.9 • 

10~ 6 

1.2 • 10~ 5 

LF-2 

-2.9- 

10^ 3 

-3.1 

10“ 3 

-2.9- 

10 -3 

1.1 ■ 10~ 3 

5.1 • 10“ 5 

-8.6 ■ 

10~ 3 

5.5 • 10~ 2 

Compressed-1 

-1.6- 

10^ 7 

-1.6 

10- 7 

-9.6- 

10^ 6 

3.4 • 10~ 2 

-1.6 • 10 -3 

-9.6 ■ 

10~ 6 

-1.4 • 10 -5 

Compressed-2 

-9.1 • 

10^ 3 

-9.4 

10" 3 

-9.1 • 

10^ 3 

3.5 • 10~ 2 

-1.5 • 10“ 3 

-2.7- 

10~ 2 

7.1 • 10~ 2 

VLF-1 

-3.6- 

10- 7 

-3.6 

10- 7 

-3.6- 

10- 7 

-3.6 ■ 10~ 7 

-4.1 • 10 _1 ° 

-1.1 • 

10~ 6 

-2.2 ■ 10 -6 

VLF-2 

-6.6- 

10~ 4 

-9.3 

10“ 4 

-6.6- 

10~ 4 

CO 

1 

0 

<N 

CN 

1 

-8.8- 10“ 5 

-2.0 • 

10~ 3 

CO 

1 

0 

0 

1 

Twin-1 

-1.9- 

10- 7 

-1.9 

10- 7 

-1.9- 

10- 7 

1.4 • 10~ 6 

-8.1 • 10“ 8 

-5.6 • 

10- 7 

-7.5 ■ 10~ 7 

Twin-2 

-0.21 

-0 

21 

-0.21 

-0.18 

6.9- 10“ 4 

-0.52 

-0.61 


TABLE V. Numerical results for the Higgs-coupling modifications for the benchmark points defined in table I. 



2HDM-3 

2HDM-4 

VLF-3 

Twin-3 

Twin-4 

A w 

-0.69 

-0.01 

-0.01 

-0.61 

-0.06 

A z 

-0.69 

-0.001 

-0.001 

-0.60 

-0.05 

A t 

1.8 

-0.05 

-0.07 

0.34 

0.06 

A b 

1.8 

-0.05 

-0.07 

0.34 

0.06 

A 7 

-0.05 

-0.09 

-0.07 

-0.03 

-0.05 


TABLE VI. Higgs-coupling modifications for the benchmark 
points of the LET model, defined in table III. 

suppressed by cos a = 0.94 and the fermion couplings 
are enhanced by 1/cosa = 1.06, and — $2 mixing is 
negligible. 

We have illustrated that the LET model is predictive 
and allows for various interesting phenomenological sce¬ 
narios. The model allows for large coupling modifications 
as well as moderate ones that can be expected to be mea¬ 
surable, hence the model is testable. A more thorough 
analysis of the parameter space is required to see which 
parameter values are preferred by experiment. 

VI. VERY LIGHT FERMIOPHOBIC SCALARS 

The VLF scenario contains a fermiophobic scalar par¬ 
ticle with a mass of O (1 GeV) in addition to the SM-like 
Higgs. Such a particle would only decay into pairs of 
photons, and is not necessarily ruled out since it could 
have escaped detection so far. We now review the exper¬ 
imental bounds that are relevant to this scenario. 

Fermiophobic Higgs particles are not unique to the 
SB model: they also appear in a type-I 2HDM with 
a = 7t/2 [44] and in models with 517(2)i-triplet Higgs 
fields [45, 46]. Mass bounds from direct searches are read¬ 
ily available in the literature [47]. Assuming SM cross 
sections, the four LEP experiments [48-51] have placed 
a lower limit m# > 107 GeV on the mass of a fermio¬ 


phobic Higgs by looking for decays into pairs of photons. 
More recent searches by ATLAS [52] and CMS [53] in 
the diplroton channel as well as the WW, ZZ channels 
[54] extend this lower limit to mn > 194 GeV. However, 
the cuts on the energy of the photon pair in the LEP 
analyses make these searches insensitive to fermiophobic 
Higgs particles with masses below 10 GeV. 

In order to see to what extent the lower mass bound 
applies to the fermiophobic Higgs 77° °f the SB model, we 
need to examine its couplings to SM particles. The field 
has no tree-level couplings to fermions, and its neutral 
component is a SM singlet. Hence H ° only couples to the 
SM through W — W' mixing, Z — Z' mixing and scalar 
mixing. The former two are negligible since they are of 
order £ 2 . Scalar mixing can become substantial in the 
Twin scenario, but it is negligible in the VLF scenario. 
Hence H% has only O (£ 2 ) couplings to the SM in the SB 
model, and therefore the experimental bounds are evaded 
trivially. 

This may change in the LET model. We have seen that 
a very light fermiophobic Higgs with a mass of a few GeV 
becomes a natural possibility. Since it is <f> 2 -like, it has 
significant couplings to W, Z , distinguishing it from the 
very light fermiophobic Higgs of the SB model, distin¬ 
guishing it from the very light fermiophobic Higgs of the 
SB model. Hence such a state could have significant pro¬ 
duction rates at the LHC. A more thorough analysis of 
the production cross-section is necessary in order to pre¬ 
dict the signal strength, which is beyond the scope of this 
work. 


VII. A TWIN HIGGS SCENARIO 

The Higgs signal strength in each channel has been 
measured at the LHC [2, 55, 56]. These measurements 
constrain the Twin scenario. In order to compare the 
Twin benchmarks to these experimental results, we have 
to consider the production cross-sections and branching 
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ratios for both ft 0 and 77°. 

We parameterize the twin states S° = h° ,77° in terms 
of the gauge eigenstates as 

hr =ai/i°n + a 2 ft° 22 + «3 h° R , 

H° 2 =a F 1 p hl n + a F 2 p h 0 h 22 + a pp ft° fl , (14) 

where FP stands for fermiophobic. The coefficients for 
the Twin-1 benchmark are given by 


(ai,a 2 ,(i 3 ) =(0.9999,0.0166,0.0006), (Twin-1) 
(af p , a 2 P , af P ) =(-0.0006, -0.00001,0.9999), (15) 

whereas for the Twin-2 benchmark they are 

(ai,a 2 ,a 3 ) =(0.7874,0.0136,-0.6162), (Twin-2) 
(af p , a 2 p , af p ) =(0.6162,0.0096,0.7876). (16) 


We need to find the signal strengths ^ X (S°) for each de¬ 
cay channel S° —> xx: 

f c0'_ vjpp s °) X BR(S° -» xx) 

^ > a{pp h o )sM x BR{h o xx)sM > y •) 

where cr{pp S°) is the production cross-section for 
S°. Then we need to add the signal strengths of ft 0 and 
77°. In order to estimate the magnitude of the deviations 
from their SM values, we neglect loop corrections in the 
following discussion. 

At the LHC, the Higgs can be produced in vector- 
boson fusion, VH associated production, gluon fusion, 
and production in association with tt pairs [57]. The 
former two processes are proportional to the Higgs cou¬ 
pling to vector bosons, whereas the latter two scale with 
the top coupling (the main contribution to gluon fu¬ 
sion comes from a top-quark loop). For ft, we have 
Ai = A w = A z = At = ai — 1 (see appendix C, we 
neglect O (£ 2 ) and O (/?i) corrections). Hence as a tree- 
level approximation we have 


a[pp —> ft 0 ) 
cr(pp h°) SM 


(18) 


Similarly, Af p = A(^ = A|) p = Af p = af p — 1 , up 
to O (£ 2 ) and O (/3i) corrections. Hence at tree level we 
have 


a(pp 77°) = F p 2 
cr(pp -» h°) S M 1 


(19) 


For the Twin-1 benchmark we have a\ = 1, (a pp ) 2 = 
4 • 10 -7 (see eq. (15)). That is, h° is produced at the 
same rate as the SM Higgs, whereas H 2 production is 
suppressed. In the Twin-2 benchmark, however, both 
states have a significant production rate since a? = 0.62, 
(a pp ) 2 = 0.38. 

As for the branching ratios, we only take into account 
the decay channels listed in table VII; all other channels 
have negligibly small branching ratios. The given SM 


values were calculated with HDecay [58] using m/, = 126 
GeV. We estimate the corresponding branching ratios of 
the SB model using the Higgs-coupling modifications. 
According to eq. (10) the b quark and the r couple to 
the same scalar gauge eigenstate, namely h { { 22 . We thus 
assume that A r = At, = A 2 and A pp = A pp = A pp . 

Now we are ready to calculate the branching ratios 
of the SB model. The partial decay widths for the ft , 0 
and H 2 decays into yiyi = WW, ZZ, gg, 77 , cc scale 
with respectively (I + A 1) 2 and (1 + Af p ) 2 , whereas they 
scale with respectively (1 + A 2 ) 2 and (1 + Af p ) 2 for the 
2 / 22/2 = bb. tt decay channels. Thus the branching ratios 
for the SB model are given by 


BR(h° —/■ xx) 

BR{H 2 -f xx) 
Blh 


(1 + A x ) 2 BR(h° -» a:x)sM 
_ (1 + Ai) 2 Hi?i + (1 + A 2 ) 2 BR 2 ’ 

(1 + A FP ) 2 BR(h° —> xx)sm 
~ (1 + A pp ) 2 Hi?i + (1 + Af p ) 2 BR 2 ’ 

= BR(h° —> 2/i2/i)smj 

yi=W,Z,s,7,c 


BR 2 = Y, BR ( h ° ->• 2 /22/ 2 )sm. (20) 

V2=b,T 


In the Twin-1 benchmark we have |Ai i2 | <C 1 (see ta¬ 
ble V), hence the branching ratios for ft 0 barely deviate 
from their SM values. The branching ratios for H 2 are 
small, since |l + Af p | = |<Ji 2 | “C 1. The branching ratios 
for the Twin-2 benchmark are listed in table VII. 

The total signal strength is given by the sum of the 
contributions from ft 0 and H 2 . Combining eqs. (18)- 
(20), we find 

O-i (1 + A ^) 2 

LLgc tot — --—— — -- -... 

(1 + Ai) 2 BRi + (1 + A 2 ) 2 BR 2 

+ (qfT(i + Ag p ) 2 _ (21) 

(1 + A fp ) 2 BRi + (1 + A FP ) 2 BR 2 ' 

We find y x ,tot = 1 for the Twin-1 benchmark: ft 0 con¬ 
tributes with the SM strength whereas H 2l despite its 
identical mass, is hardly produced and cannot show it¬ 
self by decays to SM final states. The signal strengths 
for the Twin-2 benchmark have been summarized in ta¬ 
ble VII. In this case two states with the same mass but 
different decay properties are present. Still, none of the 
decay channels has a total signal strength that deviates 
significantly from 1. The reason is that the SM contri¬ 
butions to the signal strength are simply divided among 
the two scalars. 

The situation may change in the LET model, in which 
the fermiophobic scalar 17° can be the twin partner of 
ft 0 . The state 77° is a mixture of components of both 
$1 and <b 2 . Since the latter is an SU (2 )l -antidoublet, 
it couples to W, Z. After scalar mixing, these addi¬ 
tional contributions can give the twin Higgs a total signal 
strength that differs from the SM prediction. Note that 
the Higgs-coupling modifications of the LET model are 
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ww 

ZZ 

99 

77 

cc 

bb 

TT 

BR(h° — » xx)sm 

0.216 

0.027 

0.077 

0.002 

0.026 

0.594 

0.057 

BR(h° -S- xx) 

0.206 

0.026 

0.073 

0.002 

0.025 

0.610 

0.059 

BR(H% -> xx) 

0.235 

0.029 

0.084 

0.002 

0.028 

0.566 

0.054 

/ix(h°) 

0.59 

0.59 

0.59 

0.59 

0.59 

0.64 

0.64 

MH°) 

0.41 

0.41 

0.41 

0.41 

0.41 

0.36 

0.36 

M*(fc°) + MH°) 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 


TABLE VII. Branching ratios and signal strengths of h° and 
H® decays in the Twin-2 benchmark of the SB model. The SM 
values of the branching ratios were calculated with HDecay 
[58] for rrih = 126 GeV. We neglect the branching ratios for 
the /r/i, ss, tt, Z 7 decay channels. 

not universal (see table VI). This means that the Higgs 
production cross-sections do not scale trivially, as they 
did in the SB model. Hence a prediction of the twin- 
Higgs signal strength requires a more detailed analysis of 
the production cross-section. 


VIII. CONCLUSIONS 

In this work we have studied the low-energy phe¬ 
nomenology of the trinification model as described in 
refs. [7, 19-23]. It is based on the trinification group 
SU(3)c x SU(3)l x SU(3)r. In order to simplify our 
study, we have integrated out the fields that obtain 
masses of the order of the trinification scale. This re¬ 
sulted in a left-right-symmetric model with two scalar 
bidoublets $i, $ 2 and one right-handed doublet $_r. 
While the bidoublets obtain vevs of the order of the weak 
scale, the right-handed doublet has a vev M that de¬ 
scribes the scale at which the left-right symmetry is bro¬ 
ken. It may be as low as a few TeV but could also be 
much higher. Only $1 couples to fermions. We call this 
effective model the low-energy trinification (LET) model. 
As an intermediate step towards a better understanding 
of the LET model, we have studied the Single-Bidoublet 
(SB) model, a simplified model in which $2 has been set 
to zero. 

Our Ansatz for the Yukawa sector was based on the 
Yukawa Lagrangian of the trinification model. The free 
parameters were used to fix the masses of the top and 
bottom quarks, which are the most important for com¬ 
parison to experimental searches. We left out a discus¬ 
sion of the first and second fermion generations. An im¬ 
proved version of the LET model containing these lighter 
fermions would require the introduction of new Higgs 
fields, mixings with heavy states, and the consideration 
of renormalization-group effects. 

In order to showcase the possible phenomenological 
scenarios of our model, we have defined a set of bench¬ 
mark points for both the SB model and the LET model. 
In the SB model, all new scalars decouple from the SM 


unless we tune some of the dimensionless scalar-coupling 
constants to O (v 2 /M 2 ) values. For such small cou¬ 
pling constants, interesting phenomenological scenarios 
are possible at low energies, like a fernriophobic scalar 
particle with anO(l GeV) mass in addition to a SM-like 
Higgs, higher-mass states with fernriophobic components, 
or a degenerate state (‘twin Higgs’) at 126 GeV. On the 
other hand, the full LET model always has at least one 
other CP-even scalar, one CP-odd scalar, and a pair of 
charged scalars with masses in the O (100 GeV) range. 
They appear without tuning the coupling constants to 
very small values. The aforementioned phenomenologi¬ 
cal scenarios therefore become natural possibilities. 

To show to what extent these scenarios can be dis¬ 
tinguished from the SM in experiment, we have calcu¬ 
lated the Higgs-coupling modifications for the bench¬ 
marks. For the SB model, they are negligibly small in 
most cases. However, the benchmarks of the LET model 
lead to sizable effects on the Higgs couplings. Parts of 
the parameter space can already be excluded using the 
known limits given in the literature. 

Subsequently, we have studied the scenario with a very 
light fermiophobic (VLF) Higgs in more detail. Such a 
particle decays only into pairs of photons and escapes the 
currently available bounds from direct searches. For the 
VLF Higgs of the SB model, the signal strength is very 
small since it has only O (y 2 /M 2 ) couplings to the SM. 
We argued that the VLF Higgs of the LET model may 
have significant production rates at the LHC. We also 
studied the twin Higgs scenario (degenerate Higgs state) 
in more detail. For the SB model, we expect no signifi¬ 
cant deviations from the SM prediction even though each 
of the two degenerate states has different decay proper¬ 
ties. On the other hand, the twin Higgs of the LET model 
may lead to significant deviations, because the bidoublet 
$2 can introduce direct couplings of the new state to 
SM particles. A more detailed analysis of the production 
cross-section appears necessary for a detailed comparison 
with the measured signal strengths. 

We have calculated the phenomenological scenarios of 
the LET model at several benchmark points and found 
interesting consequences for the properties of the Higgs 
bosons. Because of the number of coupling constants in 
the potential a systematic investigation of all regions of 
the parameter space is still missing. Nevertheless, the 
examples given show the large variety of possibilities still 
not excluded by experiment. The LET model turned out 
to be an interesting extension of the SM. It is predictive 
and most properties can be tested or constrained using 
forthcoming LHC data. 
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The electromagnetic coupling constant is found to be 


Appendix A: Gauge sector of the LET model 


_2 g'gLgn _ 

V 4 9 ,2 (9 2 L + c JR) + gl9R 


g L sin d w . 


(A7) 


The gauge-boson sector consists of seven fields: W RR 3 
for SU(2) Lt R and B for U(1)b-l- The fields TF|, W jj, 
and B are neutral, whereas the remaining fields mix to 
form the charge eigenstates W RR = (Wi R TiWl R )/y/2. 
These charged states are rotated by an angle ( into two 
pairs of charged mass eigenstates W ± and W r± : 


( _ / cos£ sin£)\ (w R \ 

[w^j ~ y-smc coscy \wi )' 


(Al) 


Here the W ± correspond to the charged vector bosons of 
the SM. The B- bosons are new massive vector bosons. 
The mixing angle is very small (we define s x = sin a:, 
c x = cosx for the sake of brevity): 


C — ~£ 2 (c„S(2/3i) + s 1 s (2/3 2 )) + O (C 4 ) • (A2) 

The masses of the charged mass eigenstates are given by 

m W =^2~ (y — y(4s(2/3i) + s a s (2/3 2 )) 2 + ® (^) ) i 

mw {1 + Z 2 + O (£ 4 )). (A3) 


Constraints on the W' and Z' masses as well as their 
mixing angles are widely available in the literature and 
allow us to constrain the new parameters gR , g ', and M. 
Direct searches for W' and Z' have been performed in 
various decay channels [59], but mass bounds are only 
given under the assumption that the new vector bosons 
have the same couplings to fermions as W and Z re¬ 
spectively. Hence they do not apply to the LET model. 
More general constraints come from fits to electroweak 
precision data [60, 61] and high-precision measurements 
[62-64]. The right-handed coupling gR is constrained by 
the bound gR/gL = 0.94 ± 0.09 from ref. [62], Combin¬ 
ing this bound with eq. (A7) and the experimental values 
g L = 0.65170 ± 0.00008, e = 0.313402 ± 0.000017 gives a 
constraint on g' as well: 


g R = 0.61 ±0.06, g' = 0.22 ±0.01. (A 8 ) 

The strongest constraint on M comes from the bound 
—0.00040 < 77 < 0.0026 from ref. [60] on the Z — Z' 
mixing angle. Combining it with eqs. (A5) and (A 8 ) we 
find 


The three neutral gauge fields mix into a massless photon 
A , a massive Z as in the SM, and a new massive state Z'. 
These states are obtained through a rotation over three 
mixing angles &w, d' w , 77 : 

A =S0 W Wl + Cg w (Sff> w W R + ce> w B ), 
z =c 8\ V C r] W L + (cg' w S v — Sg w SQ' w C rl )W / R 

~ { s 6w c d' w c r) + s 8' w s ri)B, 

Z 1 = — Cg w S rj Wl + ( CQ' W C V + Sg w Sg' w S v )W R 

+ { s o w co' w s v -s 0 ' w c ri )B. (A4) 

The angle dw is the Weinberg angle; 0' w is an analogon 
of dw for the breaking of the left-right symmetry; and 
?7 is the Z — Z' mixing angle. These angles are given in 
terms of the gauge couplings by 


M > 3.6 TeV. (A9) 


Appendix B: Scalar mass eigenstates 


The scalar fields (j = 1,2), can be parameter¬ 
ized in terms of gauge eigenstates as 


T, = 


= h 


( V 3 + h j, 11 + ia R 11 

V2 

h U 

M + h° R + ia° R 


'j, 12 

•j, 22 + lUi j, 22 


\ 


bj T h ^ 20 T ici® 


\/2 


°R 


V2 


(Bl) 


sin dw = 


Zg'gR 


V 4 9 r2 (g' 2 L + g 2 R ) 


9 9 

9l9r 


sin d'w — 


2 9 ' 


vW 


tan 77 = 


g 2 RV^g r2 (gl + 9 2 R)^ 

n ! 2\2 


9l9r C 2 


(9 2 r+W 2 ) 

The masses of the states Z , Z' are given by 


e+o^ 4 )- ( as ) 


After spontaneous symmetry breaking, they are mixed 
into eigenstates of the mass matrix. 


1. Single-Bidoublet model 


m z 


mz' 


9 lv 

2 cos dw 
9rM 
2 cos d' w 




£, 2 cos 4 d' w 
2 



In the SB model, only $1 and are present, con¬ 
taining twelve real scalar components in total. They are 
mixed into six Goldstone bosons and six massive parti- 
(A 6 ) cles. The scalar masses are given in terms of the model 
parameters by 



13 


(A 4 + A 5 S^j ) 2 


v h° 


V 

m H° 

M 2 


— Al + A2cjg 1 — 


A 3 

„2 ^2 \2 0 2 


o (a, 


l H° 


TO 


AO 


=A 3 + £ 2 


As 


(A 4 + A 5 S /3i ) 2 A 5 s (2/3i) C (20i) 


As 


_ ^2 

M 2 2 c (2/3l) 2 ? ’ 


m ff± 


A 5 — 2A 3 C(2 j (3 1 ) 
A 5 


M 2 20(2/3!) 

+aaY 


M 2 20(2/3!) 


(1 + ^ c ( 2 / 3 !)) • 



A 5 5 (2/3i) C (2/3i) / 2 

A 5 — 2A 3 C(2/3i) 


(B2) 


Here we define ft 0 as the scalar that is the most h® n - 
like and H 2 as the scalar that is the most /i^-like. The 
massive CP-odd state A 0 is a mixture of a® 11 and a® 22 , 
and the charged states are a mixture of hf 21 an d 
hf 12 with an C (£) admixture of /i^. 

2. Complete LET model 

If we include $2 into the scalar sector, we have 20 real 
scalar components. These are mixed into six Goldstone 


bosons, five massive CP-even states, three massive CP- 
odd states, and three pairs of massive charged scalars. 
We define the CP-even mass eigenstates h°, H °, H 2 , H ®, 
P 4 respectively as the most h { { 11 -, h 2 11 -, h\ 2 ti ^2 22 "; 
/i^j-like scalars. The massive CP-odd states A ( (, A 2 , A 3 
are dehned respectively as the most a® 2 2 - ; a 2 ,ii~> a 5 j-like 
scalars, and the charged states H 2 , Hf are defined 
respectively as the most hf 21 -, hf 12 -, h 2 l2 -\\ke scalars. 
Their masses are found to be 


=\( Aic «+ A2S “ ± 'J ( A i c “—A2S 2 ) 2 +a§s 2 2q) +o (a ) 


TO 


H° A 5 


M 2 2C(2/3i) 


+o(e) 


in 


m 


a 4 


M 2 2c 


( 2 A 0 


o (a , 


m 


Hi 


M 2 


= a 3 +o (e) 


m 


A° 


— - (A 7 + Aio)(c| 1 c| 2 + S^Sfa) ^-S(2/3i)S(2/3 2 ) + O (£“) , 


l A° As 


At 2 2C( 2 / 3i) 


+ 0 (a 




a 4 


M 2 2c ( 2/3 2 ) 


o(e a ) 


1 ~ ~ ~ A g 9 

F - = ~ 9 ( A 6 + A 7 + Aio)c 2 / 3 l _ / 3 2) - — C(2 / 3 1 )C(2/3 2 ) + O (£ 2 ) , 


4 


W 2 ± 


2 2 ' 
A 5 


o (a 




A 4 


M 2 20(2/3!) V ' ’ Al 2 20(2/3 2 ) 
Here we have defined the parameter combinations 


o (a. 


(B3) 


Ai =Ai + A2 Co — 

^3 


2 (A 4 + A 5 S ^) 2 ~ 2 ( A 3 + A 4 s | 2 ) 2 


A 2 = Ai + A2C^ 2 — 


A 3 


A 3 = - —-^^ + As + (A 6 + A 7 )c 2 / 3 i _ / 3 2) + (As + A 9 + Aio)(t^! c | 2 + SfaSpJ. 


(B4) 


Appendix C: Higgs-coupling modifications of the SB 
model 

The SB model has a larger Higgs sector than the SM. 
Since the SM-like scalar h° is a mixture of the various 


Higgs fields, its couplings generally differ from those of 
the SM. These couplings depend on the scalar mixings, 
which in turn depend on the vevs and scalar-potential 
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parameters. We give the resulting Higgs-coupling modi¬ 
fications A x , as defined in eq. (12), in the limit of small 
£. The modifications of the tree-level couplings to W, Z , 
t, b are given by 


Here A s (x) are the scalar loop functions, t x = 4 
and we defined the constant Asm = Ai(tw) + 
N c Q 2 Ai/ 2 {T t ) = —6.5. The trilinear and quartic Higgs 
self-couplings are modified as well: 


A 


W _ 2 

£2 ^ S (2/3i 


^G/h) s (4 l g 1 )S(2 < 3 1 )A 45 3 

S A r, 4 


„2 \ 2 

s /3i A 453 


+ 0(0 


A; 


*2 = 


A2S (4/3i) , S(4/3 1 )S(2/3 1 )A 4 53 


8 Ar 


„2 \ 2 

*/3i A 453 


+ 0(0 


A t 

0 


2A 2 s| i c^ 2/3i) 


Aj 


A R 


- 0(2/3,) A 453 -Tp- + o (£ 2 ) 


A b 2A 2 c% i c 2 


!(2/3l) + 2<4c ( 2 / 3 l) A 453 - ^ + O (£ 2 ) . 


(Cl) 


Here we defined A 4 5 3 = (A 4 + A 5 SIJ/A 3 . The main con¬ 
tributions to the loop-induced photon coupling of the SM 
come from the W and t loops. In the SB model, there is 
an additional contribution from the H ± loop: 


A _0^o(r ff ±)c( 2ft ) / , 1 2 a , \ 

7 "--I A 4 +A 2 (l + -s ( 2/3i) ) + A 5 c (2/3l) 


A 4 (A 4 + Asc^) 


+ 0(e) • (C2) 


= _" -'''’V _, n / y2\ 

A 3 (A 1 + A 2 c2 i )-(A 4 + A 5 S 2 i ) 2 + 

^-A 2 A3s| i c (2)9l) + (A 4 + A 5 S 2 01 ) 2 

A3(A 1 +A2c2 i )-(A 4 + A. 5S 2 i )2 )■ 

(C3) 


Note that the coupling modifications for vector bosons 
and fermions vanish in the small-£ limit, that is, the new 
physics decouples from the SM. However, the modifica¬ 
tions of the Higgs self-couplings only vanish for the region 
of parameter space where A 2 = 0 and A 4 = — A 5 s^. In 
general, the self-couplings are modified due to the more 
complicated structure of the scalar potential. 

If we look for regions of parameter space where the cou¬ 
pling modifications become substantial, it is more useful 
to write the coupling modifications in terms of the scalar 
mixing coefficients cq, as defined in eq. (14): 

A w,z =c/s 1 ai + sp 1 a 2 — 1 + O (£ 2 ) , 

A t =ai/c 01 - 1, 

A h =a 2 /s ) 3 1 - 1. (C4) 

The analogous coupling modifications A^ p of H 2 are ob¬ 
tained by substituting each cq by af p . Since j3\ = 0.0166 
is small, we can write Aw = Az = A t = A 4 and 
A= A pp = A pp , up to O (£ 2 ) and O (/3i) correc¬ 
tions. 
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